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We analyze the time reversal of waves in a turbulent medium using the parabolic Markovian 
model. We prove that the time reversal resolution can be a nonlinear function of the wavelength 
and independent of the aperture. We establish a duality relation between the turbulence-induced 
wave spread and the time-reversal resolution which can be viewed as an uncertainty inequality for 
random media. The inequality becomes an equality when the wave structure function is Gaussian. 
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Introduction 



Time reversal is the process of recording the sig- 
nal from a remote source, time-reversing and back- 
propagating it to retrofocus around the source. Time 
reversal of acoustic waves has been demonstrated 
to hold exciting technological potentials in subwave- 
length focusing, dispersion compensation, communi- 
cations, imaging, remote-sensing and target detec- 
tion in unknown environments ( see p] and references 
therein). The same should hold for the electromag- 
netic waves as well. Time reversal of electromagnetic 
waves is closely related to optical phase conjugation 
(OPC) which used to be limited to monochromatic 
waves. With the advent of experimental techniques, 
time reversal of high frequency EM waves hold diverse 
potential applications including real-time adaptive op- 
tics, laser resonators, high-power laser systems, opti- 
cal communication and information processing, im- 
age transmission, spatial and temporal filtering, spec- 
troscopy etc 0. 

Time reversal refocusing is the result of the time- 
reversal invariance of the wave equations, acoustic or 
electromagnetic, in time invariant media. The sur- 
prising and important fact is that the refocal spot 
in a richly scattering medium is typically smaller 
than that in the homogeneous medium. That is, the 
time reversal resolution is enhanced rather than ham- 
pered by the inhomogeneities of the medium. This 
sub-diffraction-limit retrofocusing is sometimes called 
superresolution and in certain regimes has been ex- 
plained mathematically by using radiative transfer 
equations. 

In the previous experimental, numerical or theoret- 
ical results the superresolution comes as a linear func- 
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tion of the wavelength but independent of the aper- 
ture. In this letter we show that in fractal media the 
resolution can be a superlinear (between linear and 
quadratic) function of the wavelength and at the same 
time independent of the aperture. The lowest achiev- 
able refocal spot size in this nonlinear regime is on the 
order of the smallest scale of the medium fluctuations. 
Below the inner scale the resolution is diffraction- 
limited while above the outer scale it is the previ- 
ously reported aperture-independent enhanced reso- 
lution [HQ. 

We will focus our analysis on the widely used 
parabolic Markovian model for waves in atmospheric 
turbulence 0. Neglecting the depolarization ef- 
fect let us write the forward propagating wave field 
E at the carrier wave number k as E{z, x) = 
^(z, x ) e l ( fcz -"t) j x £ K 2 where the complex wave am- 
plitude W satisfies the Schrondinger equation in the 
non-dimcnsionalized form 



ilk 



dz 



7 



V{z,x) °* = o 



(1) 



with Aj_ being the Laplacian in the transverse coor- 
dinates x G M 2 and V the fluctuation of the refractive 
index. Here the Fresnel number 7 equals L z k,Q X L~' 1 
with fco being the reference wavenumbcr, L z and L x 
the reference scales in the longitudinal and transverse 
directions, respectively. The notation o in eq. lfT| 
means the Stratonovich product (v.s. Ito product). 
In the Markovian model V(z, •) is assumed to be a 
<5-correlated-in-,z stationary random field such that 



(V(z,x)V(z',x'))=5(z-z') / $(0,p)e 



ip-(x-x') 



dp 



where $(k),k = (£, p) € K 3 is the power spectrum 
density of the refractive index fluctuation and, in the 
case of atmospheric turbulence, has a power-law be- 
havior in the inertial range. For simplicity of presen- 
tation we assume an isotropic power-law 



$(k) = a H \k\ 



-1 — 2H I 



|k|G(i -\O 



(2) 




FIG. 1: The time reversal process 

where Lq and £o are respectively the outer and in- 
ner scales of the turbulence and an a constant factor. 
Usually H is taken to be 1/3 in the self-similar theory 
of turbulence. We assume that the spectrum decays 
sufficiently fast for |k| ^ £q 1 while staying bounded 
for |k| < Lq 1 . 

For this model of propagation we will prove an un- 
certainty inequality where the conjugate quantities are 
the forward wave spread and the time-reversal reso- 
lution. The inequality becomes an equality when the 
wave structure function is Gaussian. This also helps 
illustrating an experimentally observed, close relation 
between the time reversal resolution and the correla- 
tion length of the scattered wave field prior to time 
reversal [ij. 



Time reversal process 

In the time reversal procedure, a source \l/o(x) 
located at z = L emits a signal with the carrier 
wavenumber k toward the time reversal mirror (TRM) 
of aperture A located at z = through a turbulent 
medium. The transmitted field is captured and time 
reversed at the TRM and then sent back toward the 
source point through the same turbulent medium, see 
Figure 1. 

The time-reversed, back-propagated wave field at 
z = L can be expressed as 

*trM (3) 

= J G(L,x, x m )G(i,x s ,x m )^ (x s )lU(x m )(2x m cix s 
i P .(x-x s )/ w(Lj p)^^) dprfXs (4) 

where Ia is the indicator function of the TRM, G the 
propagator of the Schrodinger equation and W the 
Wigner distribution function 

W(z,x,p) = J e -^G(z,x + 7 y/2,x m ) 

xG(z,x- jy/2,x m )I A (x m )dydx m . 



Here we have used the fact that time reversing of the 
signal is equivalent to the phase conjugating of its spa- 
tial component. 

The Wigner distribution W satisfies a closed form 
equation, the Wigner-Moyal equation £|, and for the 
Markovian model its moments also satisfy closed form 
equations. In particular, the mean field equation is 



d(W) 



V x (W) = Q (W) 



(5) 



dz k 

with the scattering operator Q given by 

Q/(x,p) = ^J <f(0,q)[-2/(x,p) + /(p + 7 q) 
+/(x,p- 7 q)]dq. (6) 

Eq. JSJ) is exactly solvable and the mean refocused 
field of the parabolic Markovian model can be ex- 
pressed as 

(*tr) W (7) 

= j^T2 yrfx'o?qdw*o(q)U(x')exp[iq-x] 

x exp [z(w • (x — x') — 7 Lw • q/k — 7 L|w| 2 /2fc)] 

-fc 2 /(2 7 2 ) f L>*(-s 7 w/fc)ds 
Jo 

where the structure function Z)» is given by 

£>*(x) = J $(0,q) [l-e Jx q ] dq. (8) 

Here and below / denotes Tj the Fourier transform of 
/. The main property of we need in the subsequent 
analysis is the inertial range asymptotic: 

D* (r) « C 2 y H ' , £ Q « r « L , (9) 

where the effective Holder exponent is given by 



x exp 



+1/2 for H e (0,1/2) 
1 for H G (1/2, 1] 



(10) 



and the structure parameter G* is proportional to an . 
Outside of the inertial range we have instead D» (r) ~ 



•C in and D*(r) — > D*(oo) for r 



where 



D*(oo) > is a constant. 

Let us consider a point source located at (L,Xo) by 
substituting the Dirac-delta function 5(x — xo) for 
in Q. We then obtain the point-spread function for 
time reversal Vtr — 7*oIt r with 



(|x| 2 



>I A ( — (x-x ) 



Po(x-xo) 

r 4r (x- Xo ) (ii) 

-fc 2 /(2 7 2 )L / D,(-s(x-Xo))ds . 
Jo 



exp 



2 



In the absence of random inhomogeneity the function 
T tr is unity and the resolution scale po is determined 
solely by V : 



XL 2tt 



(12) 



This is the classical (Rayleigh) resolution formula 
where the retrofocal spot size is proportional to A and 
the distance to the TRM, and inversely proportional 
to the aperture A. 



Anomalous Focal Spot-size 

First we consider the situation where there may be 
an inertial range behavior. This requires 



fc 2 7 _2 D*(oo)L > 1 



(13) 



where D*(oo) = lim r _» 00 D*(r). 

In the presence of random inhomogeneities the 
retrofocal spot size is determined by To or T tr depend- 
ing on which has a smaller support. For the power-law 
spectrum we have the inertial range asymptotic 

T tr (x) 

~ exp [-C 2 fc 2 L|x - x | 2H *7" 2 (4ii* + 2)" 1 ] 

for £q <C |x — xo| <C Lq. Under the following condition 

(akj^VI) 1 / 11 * » fry-^A ~ Po 1 (14) 

the function T tr is much more sharply localized around 
Xo than Vq. We define the turbulence- induced time- 
reversal resolution as 

Ptt = ^ J |x-X | 2 T t 2(x-X )dx/ J T t 2 .(x)dx 

which then has the inertial range asymptotic 



Ptt 



7A 



l/H, 



take the limit of vanishing Fresnel number 7 — > in 
eq. JSJ by setting x = 7y. Then we have 



limD,( 7 y) = D \y\ 

7— >o 



(16) 



D = I / 4>(0,q)|q| 2 dq. 



The resulting mean retrofocused field (^tr(y)) is 
Gaussian in the offset variable y and the refocal spot 
size on the original scale is given by 

p tr ~ 7 A(^ L)- 1 / 2 . 

Hence the linear law prevails in the sub-inertial range. 

Turbulence-induced aperture and duality 

Intuitively speaking, the turbulence-induced aper- 
ture is closely related to how a wave is spread in the 
course of propagation through the turbulent medium. 
A quantitative estimation can be given by analyzing 
the spread of wave energy. 

Let us calculate the mean energy density in 
(^(z, x)| 2 ) with the Gaussian initial wave amplitude 



*(0,x)=exp [-|x| 2 /(2a 2 )]. 



We obtain 



<l*(^x)| 2 ) 



/ a 



V/4+7 2 L 2 /(4fc 2 a 2 )] 



x exp 



-ky(2~/ z )L / D^Lws/k)ds 
Jo 



Hence the turbulence-induced broadening can be iden- 
tified as convolution with the kernel which is the in- 
verse Fourier transform T~ X T of the transfer function 



T(w) = exp 



ky(2'f)L / D*(jLws/k)ds 



to < ^ <C Lq. (15) In view of (jni), we obtain that 



The nonlinear law (|15f) is valid only down to the 
inner scale to below which the linear law prevails 
Ptt ~ 7 A£ -1 / 2 . We see that under (T5 jl -fI3 fr p tr is inde- 
pendent of the aperture, has a superlinear dependence 
on the wavelength in the inertial range. Moreover, the 
resolution is enhanced as the distance L and random 
inhomogeneities (C*) increase. This effect can be ex- 
plained by the notion of turbulence-induced aperture 
which enlarges as L and C* increase because the TRM 
is now able to capture signals initially propagating in 
the more oblique directions (see more on this below). 

To recover the linear law previously reported in [2j , 
let us consider the situation where p tY = 0(7) and 



_ 1 , , k 1 , A/^c , 

T T(x) = 

We define the turbulence-induced forward spread <r* 
as 



|x| 2 T 2 (x)dx/ / T 2 (x)dx 



which together with p tr then satisfies the uncertainty 
inequality: 



(17) 



3 



The equality holds when T tr is Gaussian, i.e. when 
Pti < A> or g <C Ptr -C io w ith ff* = 1. This strongly 
suggests the definition of the turbulence- induced aper- 
ture as A* — 7 AL / p tr which is entirely analogous to 

Because the coherence length of the wave field is 
closely related to the spread, it is not surprising then 
to find that the turbulence-induced (de)coherence 



length 5* associated with (^(L, x)\&(L, y) J is directly 

related to ptr- Indeed, one can show that 5* rs ptr 
when the effect of the turbulent medium is dominant 
over diffraction 



Discussion 

In summary, we have shown for the parabolic 
Markovian model that, first, the time reversal reso- 
lution can be aperture independent and depend on 
the wavelength in a nonlinear way. This is due to 
the self-similar nature of the media. Second, we prove 
an uncertainty inequality for random media where the 
conjugate variables are the forward wave spread and 
the time-reversal resolution. The equality is attained 
when the wave structure function T tr is Gaussian. 

The preceding analysis has been carried out for a 
narrow-band signal. Because of the linearity of the 
equation a wide-band signal Uo(t, x) can be decom- 
posed into frequency components each of which can 
be analyzed as above and then resynthesized. The 
mean retrofocused signal can be calculated as 



(w tr ) (r, x) 

1 



«o(t,y) /iUl^i^K 



(18) 

-ik(t+r) 



2tt7 2 £ 2 J ""-'"J ~ \ 7 L 

xe *fc|x|V(27i) e -*fc|y| 2 /(27i) fc 2 Tti ( x _ y)dMt 



from which it follows that the turbulence-induced 
spread in time is given by convolution with a Gaus- 
sian kernel because T tr is Gaussian in k, see 
The Gaussian kernel has a offset-x-depending variance 
of r (x) = L J Q D*(sx)ds/-f 2 which grows rapidly with 
the offset if L ^S> 1,7 -C 1. It is precisely this rapid 
change of temporal dispersion rate with the offset that 
produces the sharp spatial retrofocusing of the time- 
reversed pulse. 



Our results above have been limited to the mean 
value of the time-reversed retrofocused field. Its sec- 
ond or higher moments can be determined from those 
of the Wigner distribution which are not exactly solv- 
able. However, the mean field is sufficient for de- 
termining all the higher moments in case of self- 
averaging. Self-averaging occurs, for example, when 
the narrow-band beam width in the transverse direc- 
tions is large compared to the correlation length of 
the random medium or when the signal is wide-band 
[l| . The former case has been analyzed extensively in 
the literature (see Q| and references therein) and there 
arise several canonical radiative transfer equations as 
the self-averaging scaling limits. The latter case of 
wide-band signals has only been studied for the one- 
dimensional medium where the issue of spatial focus- 
ing does not arise, sec 3]. In the near-self-averaging 
regime the second moment of the Wigner distribution 
can be calculated and will be reported elsewhere. 
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